Mon. Not. R. Astron. Soc. OOP. [TlfTSl C2012') Printed 8 March 2013 (MN WT^ style file v2.2) 



The SWELLS survey - V. A Salpeter stellar initial mass 
function in the bulges of massive spiral galaxies 

Aaron A. Dutton^'^* Tommaso Treu^f, Brendon J. Brewer^'^, Philip J. Marshall^, 
M. W. Auger^'S, Matteo Barnabe^ David C. Koo^ Adam S. Bolton^ 
Leon V. E. Koopmans^ 

^ Max Planck Institute for Astronomy, Kdnigstuhl 17, 69117, Heidelberg, Germany 
^ Dept. of Physics, University of California, Santa Barbara, CA 931 06, USA 

'"^Department of Statistics, The University of Auckland, Private Bag 92019, Auckland 114-2, New Zealand 
"^Department of Physics, University of Oxford, Keble Road, Oxford, OXl 3RH, UK 
^Institute of Astronomy, University of Cambridge, Madingley Rd, Cambridge, CB3 OHA, UK 

^Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94035, USA 
^ UCO/Lick Observatory, Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA 
^Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA 
^ Kapteyn Astronomical Institute, University of Groningen, P.O.Box 800, 9700 AV Groningen, The Netherlands 

Accepted 2012 October 22. Received 2012 October 12; in original form 2012 June 19 

ABSTRACT 

Recent work has suggested that the stehar initial mass function (IMF) is not univer- 
sal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or 
morphological type. In this paper, we investigate variations of the IMF within indi- 
vidual galaxies. For this purpose, we use strong Icnsing and gas kinematics to measure 
independently the normalisation of the IMF of the bulge and disk components of a 
sample of 5 massive spiral galaxies with substantial bulge components taken from the 
SWELLS survey. We find that the stellar mass of the bulges are tightly constrained by 
the lensing and kinematic data. A comparison with masses based on stellar population 
synthesis models fitted to optical and near infrared photometry favors a Salpeter-like 
normalisation of the IMF. Conversely, the disk masses are less well constrained due 
to degeneracies with the dark matter halo, but are consistent with Milky Way type 
IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale 
lengths, but due to the contribution of the bulges, the galaxies are baryon dominated 
at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with 
that found for early-type galaxies of comparable stellar mass (> IO^^Mq). Our results 
suggest a non-universal IMF within the different components of spiral galaxies, adding 
to the well-known differences in stellar populations between disks and bulges. 

Key words: dark matter — galaxies: bulges — galaxies: kinematics and dynamics 
— galaxies: spiral — gravitational lensing — stars: luminosity function, mass function 



1 INTRODUCTION 

The stellar initial mass function (IMF) is a fundamental 
quantity in many areas of astrophysics. From a theoreti- 
cal standpoint, understanding the origin of the IMF from 
first principles is ess ential to develop a. comp lete theory of 
star formation (e.g., lMcKee fc Ostrikeill2007l ). From a phe- 
nomenological standpoint, the IMF is a defining property 
of any stellar population, essential for computing quantities 
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such as stellar mass from observables, and for characterising 
their evolutionary history. 

Traditionally, the main source of empirical evidence 
regarding the stellar IMF has been our own Milky Way, 
wh ere indiv i dual s tars can be identified and counted, first 
by ISalpete 3 (| 19551 ). In the decades since then, only rela- 
tively small variations of the IMF have been found within 
our own Milky Way, despite enormous variations in the 
physical conditions within star-forming regions. These re- 
sults have been generally interpreted as evidence that the 
IMF is more or less universal — i.e. it is insensitive to 
properties of the gas and dust in which stars form (e.g.. 
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iBastian et al] I2OI0I '). If indeed the IMF is universal, the 
kinematic s of spiral |Bell fc de Jond l200lh and elliptical 
:alaxies llCappellari et alJ l2006l \ as well as lensing data 



ter halo, e.g., in spired by the results of cosmological numeri- 
cal simulations jXreu fc Koopman3l2002l: Koopmans fc Treu 



Brewer et al 



2012al . hereafter SWELLS-III) rule out IMFs 



like Salpeter's implying relatively "heavy" mass-to-light ra- 
tios. in favou r of "lighter" IMFs like that measured by 
iKroupal (|200ll ) and lChabrieil (|2003h . 

However, recently a number of independent extra- 
galactic studies have found significant deviations from the 
IMF as measured in the Milky Way. The observations 
are based on a variety of in dependent techniques, rang- 
ing from gravitational lensing ("Trcu ct al."201Cf. 'Au ger et al.l 
l20ld . ISpiniello et al] [2OII, Brewer ct al. 2012a), to stel- 
lar ki ne matics of elliptical g alaxies (e.g., iDutton et al.l 
l2012al lbl: ICappellari et al.l2012al'). and to spectral diagnostics 
of stellar populations llvan Dokkum fc Conrovl I2OI0I . I2OIII : 
ISpiniello et~al]|2012l : IConrov fc van Dokkumll2012h . 

Based on these observations, it appears that the IMF 
may depend on the stellar mass of the galaxy and hence on 
the cosmological time at which the stars formed, possibly 
reflecting the evolving physical conditions in the expanding 
universe. Or perhaps it could depend on the stellar velocity 
dispersion, reflecting the depth of the potential well. It is 
also possible that the IMF might depend on the morpholog- 
ical type of the galaxy, in the sense that early-type galax- 
ies might have "heavier" IMFs than their spiral counter- 
parts. However, the dependency of t he distribution of mor- 
pholo gical types on stellar mass fe.g.. lBlanton fc Moustakaj 
l2009f l. as well as the presence of two very distinct stellar 
populations in the bulge and d isks of typical galaxies (e.g., 
IWvse. Gilmore. fc Fran^ll997l l make it unclear whether the 
more important parameter is the overall stellar mass, or stel- 
lar velocity dispersion, or the morphological type. 

In this paper we aim to gather some insight into the 
physical origin of the non-universality of the IMF by look- 
ing for variations withm individual spiral galaxies. Is the 
normalisation of the IMF the same for bulge and disk, just 
varying with total stellar mass? Or is there one universal 
normalisation for disk-like stellar populations (presumably 
Chabrier-like) and a "heavier" (Salpeter-like) one for the 
older and more metal rich stellar populations found in mas- 
sive bulges and spheroids? Clearly the two hypotheses are 
not mutually exclusive, as a combination of both could be 
at play. 

Our strategy to constrain the IMF normalisation con- 
sists of comparing stellar masses derived from dynamics and 
strong lensing to those derived from stellar population syn- 
thesis (SPS) models. The main challenge of this approach 
- and of all other dynamical approaches - is that lensing 
and dynamics are sensitive to the total mass, and there- 
fore one needs to disentangle the stellar mass from the non- 
baryonic dark matter. For disk-dominated galaxies it is well 
known that rotation curves can be fitted with a wide range 
of st ellar mass-to-light ratios, from zero to maximu m disk 
fe.g.. lvan Albada fc Sancisill 19861 : lOutton et"al]|2005l'l. the so 



called disk-halo degeneracy (see however Amorisco fc BertinI 
l2010h . For bulge-dominated systems, the situation is similar 
in principle (the so-called bulge-halo degeneracy). However, 
if the stellar mass of the bulge dominates the inner parts 
of the gravitational potential, it is possible to get an actual 
measurement (rather than an upper limit), with only mod- 
est assumptions about the density profile of the dark mat- 



20031: iTreu et al.ll20ld: lAuger et all I2OI0I : ICappellari et"ar 
2012al : see also Berlin et al.lll994 . and references therein). 

In order to study the IMF of the bulge and disk 
component we apply the lensing and dynamical technique 
to the sample of five massive (V2.2 ~ 250 — 30Qkms~^, 
M^^t lO^^Mp)) spir al lens galaxies in the SWELLS 
Survey (jTreu et al.il201ll . hereafter SWELLS-I) with signif- 
icant bulges (bulge-to-disk ratio > 0.5 assuming a universal 
IMF) and for which ionised gas kinematic data are available. 
SWELLS is a dedicated survey to the study of lensing spiral 
disk galaxies. Its main properties are summari zed in Sec- 
tion 2 . More details can be found in the paper by I Treu et al.l 
(|201ll ). Thus, our mass models are constrained by the pro- 
jected mass within the Einstein radius measured by strong 
gravitational lensing, by the rotation curve of the outer disk 
~ 1 — 3 scale lengths from ionised gas kinematics, and the 
high resolution surface brightness profiles which are decom- 
posed into a de Vaucouleurs bulge and an exponential disk. 
We derive stellar masses of the bulge and disk for a vari- 
ety of assumptions about the structure of the dark matter 
halo, and compare the results with the estimates from stel- 
lar population synthesis models to infer the normalisation 
of the IMF independently for the two components. 

This paper is organised as follows. In §2 we present the 
observational constraints. In §3 we present the mass models. 
In §4 we present of our main results. In §5 we discuss possible 
sources of systematic errors. In §6 we conclude with a brief 
summary. Throughout, we assume a flat ACDM cosmology 
with present day matter density, Q^n ~ 0.3, and Hubble 
parameter. Ho = 70kms~^ Mpc~^. 



2 SAMPLE SELECTION AND DATA 

In this section we describe the sample selection and give a 
brief description of the data used in this pap er. The reader 
is refer red to papers SWELLS-I, SWELLS-II (|Dutton et all 
l2011bl ). SWELLS-III, and SWELLS- VI (Button et al., in 
preparation) for more details on the SWELLS selection and 
data. 



2.1 Sample selection 

As detailed in paper I, the parent SWELLS sample is se- 
lected from the Sloan Digital Sky Survey (SDSS) spectro- 
scopic database based on the detection of emission lines at 
multiple redshifts within the 3" flber and on having photo- 
metric axis ratio < 0.6. Multiband Hubble Space Telescope 
imaging is used to confirm the lensing hypothesis. The sam- 
ple of twenty confirmed lenses spans a broad range of bulge 
to total stellar mass ratio (from 0.1 to 0.9) for a fixed uni- 
versal IMF, and a range in stellar mass of over one decade 
(~ lO^Me to IO^ '^Mq; for a Chabrier IMF). 

An important question is whether the SWELLS sample 
is representative of the overall population of spiral galaxies 
within the same ellipticity and stellar mass limits. Paper 
I shows that the size-stellar mass relation of the SWELLS 
sample is indistinguishable from that of a control sample 
of SDSS-selected galaxies using the same criteria. Paper VI 
will investigate the selection function in more detail using 
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Figure 1. Postage stamp images of the galaxies analysed in this paper. The images are obtained from multieolour HST images as 
described in SWELLS-I and III. 



Table 1. Basic properties of galaxies analyzed in this paper. All errors correspond to la. Col. 1 lists the lens ID; Col. 2 the redshift of 
the deflector; Col. 3 gives the velocity dispersion of the Singular Isothermal Ellipsoid (SIE) lens model from paper III; Cols. 4 and 5 give 
half-light sizes for the bulge (circularized) and disk (major axis) from papers I and III; Cols. 6, 7 and 8 give the stellar masses for the 
bulge, disk, and total, assuming a Chabrier (2003) IMF; Cols. 9 and 10 give the model (dust reddened) rest-frame V-band luminosities 
for the bulge and disk; Cols. 11 and 12 give the model (B-V) colors for the bulge and disk. 



ID za cTsiE flbuigo -Rdiak log{Mspl!'b) l°g(A/g'|;=d) i°fi(M§ri'^t) log{iv,b) log(iv,d) (S - ^)b (B - V)a 

[kms-l] [kpc] [kpc] [Mq] [Mq] [Mq] [Lq_v-1 [Lq.v] 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


J0841+3824 


0.116 


251.2 ± 4.4 


2.42 


23.43 


11.05 ± 0.09 


11.23 ± 0.09 


11.45 ± 0.07 


10.54 ± 0.04 


10.84 ± 0.04 


0.79 


0.70 


J0915+4211 


0.078 


195.7 ± 2.2 


1.56 


4.15 


10.60 ± 0.09 


10.17 ± 0.09 


10.74 ± 0.08 


10.03 ± 0.05 


9.73 ± 0.05 


0.81 


0.72 


J0955+0101 


0.111 


238.4 ± 7.3 


1.62 


3.72 


10.63 ± 0.09 


10.17 ± 0.09 


10.76 ± 0.07 


10.00 ± 0.04 


9.73 ± 0.04 


0.87 


0.73 


J1251-0208 


0.224 


203.0 ± 2.6 


1.68 


12.01 


10.68 ± 0.07 


10.96 ± 0.07 


11.14 ± 0.06 


10.16 ± 0.04 


10.68 ± 0.03 


0.82 


0.65 


J1331+3638 


0.113 


248.1 ± 4.4 


2.86 


12.46 


10.89 ± 0.10 


10.46 ± 0.10 


11.03 ± 0.07 


10.38 ± 0.04 


10.34 ± 0.04 


0.78 


0.47 



kinematic data to construct the Tully-Fisher, Faber Jackson 
and Fundamental Plane correlations. 

The subsample analysed in this paper is selected from 
the SWELLS survey to include all the massive spiral galax- 
ies with significant bulges for which gas-kinematic rotation 
curves are available. We require the stellar mass in the disk 
to be less than two times the stellar mass in the bulge, a 
criterion that is met by all but two of the lenses in our sam- 
ple (we also exclude one system that has been shown to 
have a pseudo-bulge) although only 5 of these systems also 
have gas rotation curves available. A montage of our selected 
subsample of galaxies is shown in Figure [1] The sub-sample 
includes some (but not all) of the most massive SWELLS 
lenses. Quantitatively, the subsample spans a range in lens- 
ing velocity dispersion of 196-251 kms~^ , i.e. galaxies that 
are more massive than our own Milky Way, and spanning the 
230kms~^ threshold below which paper III concluded that 
a global Salpeter IMF is on average ruled out by the data. 
Basic properties of the subsample analysed in this paper are 
given in Table [T] 



Here Da is the angular diameter distance from the observer 
to the source, Dd is the angular diameter distance from the 
observer to the deflector (i.e., the lens), and Dds is the angu- 
lar diameter distance from the deflector to the source. The 
critical density thus depends only on the distances to the 
lens and source, which are known for all our lens systems 
from the SDSS redshifts and our adopted cosmology. 

For comparison to kinematics it is convenient to express 
the lensing results in terms of circular velocity or velocity 
dispersion, where Vsie = V^asm- However, since the mass 
proflles are in general not isothermal, the conversion from 
projected mass into circular velocity at the Einstein radius is 
non trivial. For a total mass profile steeper than isothermal, 
which is the case for our galaxies, the circular velocity at the 
Einstein radius is larger than the nominal circular velocity 
of the singular isothermal ellipsoid (SIE) model. 

Thus in our mass models in this paper we fit directly 
for b, rather than the mass projected within b (due to co- 
variance between Mproj(6) and fe), or the derived circular 
velocity /velocity dispersion at the Einstein radius. 



2.2 Strong Lensing 

Strong lensing models for our sample galaxies are presented 
in SWELLS-HI. The parameter that is most directly con- 
strained b;j;_strong lensing is the circularised Einstein radius 
b (|Treull2O10l . and references therein). This is deflned to be 
the radius inside which the average surface density is equal 
to the critical density for strong lensing: 



2.3 Gas kinematics 

We obtained major axis long-slit spectra for all five galax- 
ies with the Low Resolution Imaging Spectrograph (LRIS; 
lOke et all Il995l ) on the Keck I 10-m telescope between 
November 2008 and April 2011. Typical exposure times were 
60 minutes, with seeing conditions of FWHM ~ 1.0 arcsec. 
On the red side we used the 600/7500 line grating which 
gives a pixel scale of 1.26 A px~^. With a l"width slit the 
resulting spectral resolution is ~ 4.2 A, corresponding to ve- 
locity dispersion resolution of ~ 70kms~^ at the wavelength 
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of Ha. We adopt spatial samplings of ~ I'.'l — 1'.'5 (5-7 pixels) 
corresponding to ~ 1 data point per seeing FWHM. 

We measured rotation curves by fitting Gaussians to 
one-dimensional spectra extracted along the slit. To improve 
centroiding accuracy we fitted the redshift of neighbouring 
emission lines simultaneously. The primary set of emission 
lines are Ha 6563, [NH] 6583, 6548. The line flux ratios be- 
tween the two [NH] lines were fixed at the expected values. 
The velocity dispersions of the line pairs were also imposed 
to be the same, but were allowed to be different than Ha. 
The continuum was fitted with a 2nd order polynomial. 

The observed rotation velocities are lowered from the 
true circular velocities by several effects. Firstly, there is 
a sin(j) term due to the inclination, i, of the galaxy with 
respect to the line of sight. Secondly, there is the (poten- 
tial) position angle offset of the slit from the kinematic ma- 
jor axis. Thirdly there is beam smearing due to finite slit 
width and seeing. These effects are taken into account in our 
modelling. We create model velocity fields and extract flux 
weighted rotation velocities inside a i"width slit, which has a 
position angle offset from the major axis. For this calculation 
we assume the line emitting gas traces the stellar disk. This 
assumption is the major source of systematic uncertainty in 
our models. To minimise the impact of the uncertainties in 
the beam smearing model we exclude from our fits the in- 
ner few arcseconds of the kinematic data. We also exclude 
data points where the rotation curve is strongly asymmet- 
ric (such as due to the bar in J084H-3824), or where the 
signal to noise is low (such as in inter spiral arm regions). 
We treat the disk inclination angle and slit position angle 
offsets as nuisance parameters, allowing them to vary over 
suitably chosen small ranges. We determine an initial guess 
for the inclination from the disk minor-to-major axis ratio 
gd as suming cos(0 = V(9d " 9o)/(l - ^o)- where go ^ 0.2 
fe.g.. lHall et alll2012l ) is the intrinsic disk thickness. 

2.4 Bulge and disk structural parameters 

The light profiles are decomposed into bulges and disks as 
described in papers SWELLS-I and HI. The bulge light pro- 
file is assumed to follow a Sersic ti = 4 (de Vaucouleurs) 
profile, while the disk profile is assumed to be Sersic n = f 
(exponential) profile. We note that 084H-3824 is a barred 
spiral galaxy. However, the bar is significantly fainter than 
the bulge (« 2 magnitudes) and therefore it can be con- 
sidered negligible from a dynamical standpoint, within our 
desired level of precision, thus simplifying significantly the 
analysis. 

2.5 Stellar Population Synthesis Masses 

The photometry (typically in 4 bands, BVIK ) is used to 
compute SPS masse s using the method of [Auger et al.l 
(|2009l ) together with iBruzual fc CharlotI (|2003D SPS mod- 
els. The S PS masses for th e bu lge and d i sk are computed 
using both lChabrieil (|2003D and ISalpeteJ (|l955h IMFs. We 
consider stellar populations described by 5 parameters: the 
total stellar mass Msps, the population age A, the exponen- 
tial star formation burst timescale r, the metallicity Z and 
the reddening due to dust, tv- We employ a uniform prior 
requiring 9 ^ logio(Msps/MQ) < f3, the age is constrained 



such that star formation began at some (uniformly likely) 
time between f ^ z ^ 5, r has an exponential prior with 
characteristic scale f Gyr, and we impose uniform priors on 
the logarithms of the metallicity and dust extinction such 
that —4 ^ logio^ 5; —1.3 and —2 ^ logiorv ^ 0.3. We note 
that the priors are the same for the bulge and the disk com- 
ponents but are sufficiently conservative that they do not 
bi as our results . The p osterior PDF is sampled as described 
in lAuger et all (|2009l V The stellar mass estimates are the 
mean and standard deviation of the marginalised posterior 
PDF. 

As discussed by several authors, especially for the old 
and mostly dust free populations found in massive bulges, 
the uncertainties in stellar mass inferred from colors are 
dominate d by the IMF normal ization (see, e.g., recent dis- 
cussion bv lNewman et al]|2012l . and references therein). The 
degeneracies between ages and metallicities typical of color- 
based i nferences mostly can c el out when infer ring stellar 
masses (|BeU fc de JonjbOOll : [Auger et al.ll2009l '). Similarly, 
the effects of using different stellar population synthesis 
models are typically of order 0.05-0.1 dex for stellar pop- 
ulatio ns not do minated by thermally pulsating AGB stars 
(e.g.. lTreu et al.ii2dIol ). 



3 REFERENCE MASS MODEL 

In this section we define our reference model, which we will 
use in Section [4] to interpret the data. This is a very gen- 
eral model and should be regarded as the one providing the 
definitive results of this paper. However, in order to test the 
robustness of our conclusions, in Section[5]we define and ap- 
ply to our data a full battery of alternative models, includ- 
ing ones with adiabatic contraction and expansion, as well 
as models characterised by dark matter cores, by different 
implementations of the cosmologically motivated dark mat- 
ter halos introduced here, and models without dark matter. 
As we will show, our conclusions are robust with respect to 
the choice of model. Therefore, the reader pressed for time 
can skip Section [S] and jump directly to the Summary after 
Section [l] 

3.1 Overview 

Our reference model is one with two baryonic components, 
a spherical bulge and a thin disk, and a generic dark matter 
halo, that includes as a subset standard profiles motivated 
by cosmological simulations. We refer to this model as the 

"free" dark matter model. 

We model the bulge with a iHernauistI (|l990l ) profile, 
which is parametrised by its mass, Mbuigc, and a half-mass 
radius, rbuigc. We model the disk with an exponential profile, 
which is parametrised by its mass, Mdisk, and a half-mass ra- 
dius rdiak . We model the dark ma t ter ha lo with a generalised 
spherical iNavarro. Frenk. White! (|l997l . NFW) profile: 

p(r)oc (r/re)-^(l + r/rs)-^+^ (2) 

where 7 is the inner slope (NFW corresponds to 7 = 1), and 
Ts is the scale radius. The normalisation is determined by 
the virial velocity, V200, which is defined at a radius, -R200, 
enclosing a mean density of 200 times the critical density 
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of the Universe at redshift zero. The relation between virial 
velocity, virial radius, and virial mass is thus given by 

1/3 



V200 



R2 



fkms" 



(3) 



' [ft-ikpc] V [h-^MQ 

where h = Ho/WO kms'^Mpc"^ and G ^ 4.302 x 
10~® km'^s"^ kpcM©"'^. This reference model thus has 7 pa- 
rameters (4 for the baryons, 3 for the dark matter). 

3.2 Constraints and priors 

3.2.1 Baryons 

The disk and bulge sizes are fixed to the values ob- 
tained from our photometric bulge-disk decompositions 
(from SWELLS I and III). This assumes negligible mass- 
to-light ratio gradients in the bulge and disk and leaves 2 
free parameters, the masses of the bulge and the disk. For 
these we adopt uniform priors in logiQ(Af) over an interval 
bracketing the range of plausible values as inferred from SPS 
masses. In practice, we adopt a lower limit of half the SPS 
mass assuming a Chabrier IMF, and an upper limit of twice 
the SPS mass assuming a Salpeter IMF. For old stellar pop- 
ulations, this upper limit corresponds to a power-law IMF 
with slope ~ —3. 



3.2.2 Dark matter 

For the inner slope of the dark matter density profile we 
adopt a uniform prior in the interval ^5 7 < 2. The lower 
limit corresponds to a cored halo, while the upper limit cor- 
responds to isothermal, which mimics strong halo contrac- 
tion. 

For the scale radius of the halo we adopt a prior based 
on the concentratio n mass relation of ACDM haloes from 
iMaccio et al.l (|2008l ). The median relation is given by 



logio(c) = 0.830 - 0.098 log 



M200 



1012/i-iMp 



(4) 



The concentration, c, is defined to be the ratio between the 
virial radius and the radius where the slope of the density 
profile is —2: c = i?2oo/7'-2. For an NFW halo r_2 — r^, but 
for a generalised NFW halo r_2 = '"s/(2 — 7). We adopt a 
Gaussian prior on logiQ(c|M2oo) with standard deviation of 
0.11 dex, which is the scatt er in halo concentr ations found 
in cosmological simulations (|Macci6 et al ]|2008l ). This choice 
of prior is not critical to our results. Its main purpose is 
to introduce more freedom in the dark matter model than 
obtained by assuming a fixed scale radius (e.g.. FiVeu et al.l 
I2OI0I : [Cappellari et al.ll2012al ). and to ensure the scale radii 
are not unphysically large or small. 

For the virial velocity of the dark matter halo we adopt 
a uniform prior in the interval: ¥200^ ^ V200 ^ V2oo''kms~^. 
The lower limit to the virial velocity is obtained by assuming 
the baryon fraction (inside the virial radius) is equal to the 
cosmic baryon fraction, /bar = 0.16, thus 



^200 — 



GA/total 



/bar/ (1 — /bar) 



^/■^ / lift /I \ / ][,■[ . \ 



116.4 

km s^i 



10"M 



(5) 

where the total stellar mass Aftotai = Afbuige + Afdisk. Simi- 
larly, the upper limit VS)o^, corresponds to a baryon-to-stars 



conversion efficiency of 1%. This conservative limit is moti- 
vated by satellite kinematics, weak lensing and halo abun- 
dance matching, wh ich find conversion e fficiencies of > 2.5% 
in massive galaxies (|Dutton et al.|[2010l ). 



3.3 Mass model fitting 

The model, parametrized by 9, is fitted to the kinematic 
and lensing data simultaneously using a Bayesian Markov 
Chain Monte Carlo approach similar to that used in previous 
SWELLS papers. Specifically we fit for the lensing Einstein 
radius, b, and the observed gas rotation curve Vrot{ri) using 
the standard likelihood function: 



£{e) = exp 



[b ~ 6-°d(e)]^_ _ ^ [Kot(rO - K!sr^{e)r 



2al 



(6) 

As described above, our model for the gas rotation curve 
takes into account beam smearing due to seeing and finite 
slit width. 



4 RESULTS 

The results presented in this section are based on models 
with the "free" dark matter halo described in Section 3. As 
discussed in § (5] none of our main results are sensitive to 
the functional form of the dark matter halo that we adopt. 

The results of our fits to the rotation curves and lensing 
Einstein radii are shown in the right panels of Fig. (2] The 
black shaded region shows samples of acceptable models, 
which should be compared to the red and blue data points 
(which correspond to the receding and approaching sides of 
the rotation curve, respectively). The inset panels show the 
1 and 2(7 constraints on the Einstein radius from strong lens- 
ing (SWELLS-III). The model Einstein radii are shown as 
histograms. The orange shaded regions show the 68% con- 
fidence regions on the total circular velocity. The magenta, 
blue and grey shaded regions show the decomposition into 
bulge, disk, and dark matter, respectively. Overall the mod- 
els fit the data remarkably well and allow us to answer some 
specific questions. 



4.1 Is the IMF universal witiiin individual 
galaxies? 

The posterior probability distribution functions (PDFs) for 
the bulge and disk masses are shown in the left panels 
of Fig. [2] next to the corresponding rotation curves. The 
marginalised constraints on the bulge, disk, and total stellar 
masses are given in Table [5] (for reference, the correspond- 
ing masses from SPS models are given in Table [ij . In the 
reference "free" models, the disk stellar masses are some- 
what loosely co nstrained due to the well-kn own disk-halo 
degeneracy (e.g. . Ivan Albada fc Sancisill 19861 ). However, the 
bulge masses are well constrained because the inner part of 
the mass density profile is too steep to be described by the 
dark matter halo. 

The grey circles in Fig. [2] correspond to a universal 
Chabrier IMF. Three of the five galaxies in our sample 
are inconsistent (at greater than 2a level) with a universal 
Chabrier IMF in both the bulge and disk. A heavier IMF is 
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Figure 2. Mass model results for a free dark matter halo. Left: Joint constraints on bulge and disk masses. The contours enclose 68% 
and 95% of the probability. The red and blue histograms show the marginalised probability distributions for the bulge and disk masses, 
respectively. The dashed and dotted lines correspond to the median SPS masses assuming Chabrier and Salpeter IMFs, respectively. 
For reference, the grey ellipse shows the Itr uncertainties on these SPS masses assuming a universal Chabrier IMF. Right: Data model 
comparison. The black shaded region shows the 68% confidence region of the "observed" model (including inclination and beam smearing 
effects). This should be compared to the red and blue points, which correspond to the observed rotation velocities from the receding 
and approaching sides of the galaxy, respectively. In the inset panels histograms show the distribution of model Einstein radii, while the 
shaded bands correspond to the 1 and 2cr constraints on the observed Einstein radii. These grey bands are repeated in the main panel. 
For reference the black point shows the circular velocity of the SIE model (note: this is not used in the fit). In the main panels, the 
orange shaded region shows the 68% constraint on the total circular velocity profile. The magenta, cyan, and grey shaded regions show 
the circular velocity profiles due to the bulge, disk, and dark matter halo, respectively. The bulge masses are well constrained due to the 
declining nature of the circular velocity curves, while the disk masses are less well constrained, due to degeneracies with the dark matter 
halo. 
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Table 2. Summary of stellar masses from fits to lensing and dynamics data with a free dark matter halo, and a comparison with SPS 
based masses. Cols 2,3, and 4 give the stellar masses (median together with 68% posterior probability) of the bulge, disk and total 
(bulge+disk); Cols 5,6, and 7 give the IMF mismatch parameter a = Mm/Afg'pg'' (mean and standard deviation); Col 8 gives an upper 
limit to a within the Einstein radius using results directly from SWELLS-III. 



Name 
(1) 


log(J\/LD, bulge) 
(2) 


log(MLD,disk) 

(3) 


log(AfLD,total) 

(4) 


log(abulgo) 
(5) 


log(adisk) 
(6) 


log(cJtotal) 

(7) 


log(aicns) 

(8) 


J0841+3824 
J0915+4211 
J0955+0101 
J1251-0208 
J1331+3638 
Ensemble 


q-,+ 0.02 
li.O(„Q Q3 

iU./Ci_Q Qg 

io.9oin^ 

ii.29ini 


ll.UD_Q (jg 

10.29«i? 
11.16^0-? 


11 54+0.05 
in 04+0.03 

11 14+0.03 

.,7+0.05 

-0.04 


0.31 ±0.10 
0.13 ±0.11 
0.41 ±0.12 
0.18 ±0.15 
0.40 ±0.11 
0.29 ±0.05 


-0.15 ±0.13 
0.27 ±0.24 
0.12 ±0.26 
0.15 ±0.17 
0.06 ±0.25 

-0.01 ±0.12 


0.10 ±0.08 
0.18 ±0.09 
0.37 ±0.09 
0.17 ± 0.14 
0.34 ± 0.09 
0.24 ± 0.04 


0.29 ±0.10 
0.09 ±0.10 
0.41 ±0.08 
0.33 ±0.07 
0.42 ±0.09 
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Figure 3. Comparison between bulge (red, filled circles), disk 
(blue, open circles) and total (black, squares) stellar masses from 
lensing plus dynamics, Mldj with those from stellar population 
synthesis models assuming a Chabrier IMF, Mgps. The error bars 
enclose 68% of the posterior probability. The disk masses of our 
five galaxies are consistent with a Chabrier IMF, but the bulge 
(and total) masses favour an IMF a factor of ~ 2 heavier than 
Chabrier. 



required for at least the bulge. Since Chabrier IMF seems to 
be in general preferred for starforming galaxy disks, our data 
indicate that the IMF may be non-universal within galaxies. 

A comparison between the bulge and disk masses from 
our lensing plus dynamics analysis, Mld, with those from 
SPS models, Msps, is shown in Fig. [3] For the disks, 
the LD masses are consistent with a Chabrier IMF, in 
agreement with prev i ous studies of galact i c disks (e.g 
Beh fc de Jon3 I2001I: 

However for the bulges, the 
~ 2 higher than predicted 
and consistent with a Salpeter-like 
latter result is in agr eement with 



2011al : iBarnabe et all I2OI2I '). 

LD masses are a factor of 
by a Chabrier IMF 

normalisation. This ^ 

studies of massive early- type galaxies (e.g., lAuger et al.l 



interesting to note that the surface densities of the SWELLS 
bulges are comparable to the highest density ETGs in the 
local universe, which also favour I MFs with Salpeter-type 
normalisation l|Dutton et al.|[2012al ). 

Fig, m shows the P DFs for the "IMF mismatch param- 
eter" iTreu et al.] (|2010l ') a = Mld /Msps quantifying the re- 
lation between between LD masses and SPS masses. Uncer- 
tainties on both measurements are taken into account when 
deriving the posterior distribution of a. The mean and stan- 
dard deviation of the PDFs from Fig. |4] are given in Tabled 
On average, the IMF of bulges is a factor of ~ 2 heavier than 
Chabrier (logabuigo = 0.29 ± 0.05), while the IMF of the 
disks is consistent with Chabrier (logodisk ~ —0.01 ± 0.12). 
We note that in 4 out of 5 galaxies the disk masses are also 
consistent with a Salpeter IMF. The total (bulge + disk) 
masses have logototai = 0.24±0.04, which is consistent with 
a Salpeter-like normalisation of the IMF. 

In SWELLS-III we concluded that a global IMF with 
Q above the Salpeter value is ruled out (at 98% confidence) 
for galaxies with lensing velocity dispersions below 230 km 
s~^. This might seem at odds with the conclusions of this 
paper. However, 3 out of the 5 galaxies studied here have ve- 
locity dispersion greater than 230 km s~^ (see TabIe[T]). The 
galaxies with higher velocity dispersions all favor a bulge 
Q values higher than the Salpeter value, while the galaxies 
with lower velocity dispersions favor bulge a values lower 
than the Salpeter value. 

In addition. Table [2] gives the a parameter inside 
the Einstein radius calculated using results presented in 
SWELLS-III. Since these a axe based on total masses, they 
are upper limits. Thus for consistency they should be larger 
than the a we derive in this paper. A direct comparison is 
complicated because the a's are measured in different aper- 
tures, but nevertheless, the a values for the bulges that we 
derive here are fully consistent with the upper limits from 
SWELLS-III. 



4.2 Are disks sub-maximal? 

The rotation curves of spiral galaxies have contributions 
from the stellar mass, gas mass, and dark matter. The 
relative contributions of each are difficult to estimate 
from rotation curve modelling due to well-known degen- 



© 2012 RAS, MNRAS 000.[TlfT5l 



The IMF of massive bulges 9 



Bulge 
Disk 



10 




10 



O 



-0.3 0.3 

logioO( = logi(,(Mj,o/Msps) 



0.6 



5 - 



1 


1 ' 
Total 


1 


1 ' 


1 




- free 














U I 












Q) ! 






q 




- 


^ \ 






n 


- 










1 








1 




II 




- 








X 


- 




841 
















1 955 






























— 1 — 


_ . . 






1 





-0.3 0.3 

logioa = logio(MLD/Mgps) 



0.6 



Figure 4. Posterior distributions for the IMF mismatch parameter relative to a Chabrier IMF, a. Left: for bulges (red, solid lines) 
and disks (blue, long dashed lines); Right: for total (bulgc+disk) stellar mass. The bold histograms show the joint constraints on the 
mean a. Disks are consistent with a Chabrier IMF, while bulges require an IMF roughly twice as heavy as a Chabrier IMF. The overall 
normalisation is close to Salpeter. 



eracies. This led to the m aximum disk hypothesis (e.g., 
Ivan Albada fc Sancisil Il986l ) which proposes that the stel- 
lar disk makes the maximum possible contribution to the 
rotation curve. 

According to the standard definition (|Sackettlll997l ). a 
disk is maximal if its contribution to the circular velocity 
at 2.2 disk scale lengths is 0.75 < Vdisk/V2.2 < 0.95. This 
is shown as a blue shaded region in Fig. [5] The blue circles 
in Fig. [S] show that our five galaxies have Vdisk/V2.2 < 0.7, 
and are thus sub-maximal. This is in broad agreement with 
results from the Disk Mass Project (solid and dashed lines in 
Fig. [5]) which argued t hat all galaxy disks are sub-maximal 
llBershadv et al. l 2011^. in addition t o earlier studies (e.g., 
lBottemalll993l : ICourteau fc Rixlll999l) . 

The galaxies in the Disk Mass Project are disk- 
dominated, so sub-maximal disks imply their galaxies are 
not baryon dominated inside 2.2 disk scale lengths (i.e., 
there is significant dark matter). However, in our galaxies 
there is a substantial bulge contribution (usually more than 
the disk) to V2.2. The contribution of the stars (bulge plus 
disk) to V2.2 is — 0.75 to 0.95 (black points in Fig. [5]l, and 
thus these galaxies can be considered maximal in their to- 
tal baryonic content at 2.2 disk scale lengths. (See also the 
results of the analysis of a disk galaxy lens including stella r 
kinematic constraints in SWELLS-IV. lBarnabe et al.ll20l3 .) 



4.3 Are dark matter haloes cuspy? 

The use of general dark matter halos allows us to inves- 
tigate the inner slope of the dark matter density profile, 
which is predicted by numerical simulations to be approx- 
imately unity, in the absence of baryons. Fig. [6] shows 
the PDFs of 7 as inferred from our data. For individ- 




0.2 
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Figure 5. Contribution of disk, bulge and bulge-|-disk to circular 
velocity at 2.2 disk scale lengths. Di sks are sub- maximal i n agree- 
ment with the Disk Mass Project llBershady et al.ll201ll) . where 
the solid and dashed lines correspond to their best fit relation and 
scatter, respectively. However, the contribution of the disk plus 
bulge at 2.2 disk scale lengths is maximal, and thus our galaxies 
are baryon dominated within 2.2 disk scale lengths. 



ual galaxies 7 is only weakly constrained, reflecting the 
fact that the bulge stellar mass dominates in the inner 
regions. Taken as an ensemble, and assuming a univer- 
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Figure 6. Posterior distributions for tlie central slope of the dark matter density profile, 7 (left panel), and halo response parameter, v 
(right panel). The bold histograms show the joint constraint on the mean. Haloes with 7 > 1 or halo contraction u > are, on average, 
disfavoured. 



sal value for 7, values larger than unity (i.e., NFW) are 
marginally disfavoured [p(7 > 1)=0.07], consistent with 
what is found for early-type g alaxies of corn parable mass 
jTreu fc Koopmans 2004 : Duttonet al.l|2012bl), albeit ther e 



are counter- examples (fSonnenfeld et al.ll2012l : lGrilldl2012l ). 
which may indicate that there is a broad scatter stemming 
perhaps from different formation histories. 



At face value this result would imply a marginal con- 
flict with the universal profiles predicted by cold dark mat- 
ter only numerical simulations. However, since the central 
regions of these galaxies are baryon dominated, an accu- 
rate comparison depends crucially on how baryonic effects 
alter the underlying dark matter halos. Standard argu- 
ments suggest that baryonic cooling to form stars should 
steepen the overall mass density p rofile, thus causing the 
dark matter halo to steepen as well (|Blumenthal et al]| 19861 : 
iGnedin et al]|2004l ). This would only exacerbate the ten- 
sion between our observations and theoretical prediction. 
However, in practice other processes may occur that act 
to expand the halo, such as mass outflows due to stellar 

Gilmord [2OO5I : 



Read 



and/or ac tive galactic nu clei (e.g., 

iDuffv et al. 2010; Pontzcn fc Governatdl2012l )" or dynamical 
fricti on between infalling galaxies and the dark matter halo 
(e.g.. lEl-Zant et al.|[200ll : [Johansson et al.ll2009l '). Investigat- 
ing this complex physics in detail goes beyond the scope of 
this paper, but it is interesting to notice that our data sug- 
gest that the net effect of baryonic physics appears to lead 
to real halos that are flatter or at most as flat as NFW, 
and not as cuspy as the standard contraction recipes would 
imply. 



Table 3. Summary of priors on dark matter halo for our various 
mass models. U{a, b) is a uniform prior with limits a, and b. 5(a) 
is a delta function prior at a. 
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5 TESTING SYSTEMATIC EFFECTS 

In this section we discuss how our results depend on the 
functional form of the dark matter halo (ij IS.lf) . and we test 
potential systematic uncertainties arising from contamina- 
tion of the lensing signal by mass along the line of sight 
(§ I5.2|l . In § 15.31 we discuss briefly the effects of neglecting 
cold gas in our lensing and dynamical analysis. We find that 
our conclusions are robust with respect to the choice of mass 
models and line of sight effects, and the inclusion of cold gas. 



5.1 Alternative dark matter halo models 

Our inferred dark matter profile is less cuspy than the NFW 
profile and one may be concerned that we have fitted away 
part of the central dark matter into our bulge component, 
thus leading to the preference for Salpeter-like IMFs over 
Chabrier. We therefore also consider a wide range of less fiex- 
ible dark matter haloes, ranging from models with no dark 
matter or "mass-follows-light" ("mfl"), to cored dark matter 
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Figure 7. Posterior distributions for the IMF mismatch parameter relative to a Chabrier IMF, a, for eight different assumptions about 
the functional form of the dar k matter halo: mfl — mass-foUows-light (i.e., no dark matter halo); core — c ored dark mat t er hal o; nfw — 
iNavarro. Frenk. Whit'3 dark matter halo, constrained to follow the concentration-mass relation from lMaccio et ah] |200 8|) a l lowin g 

for scatter; ac — adiabatically contracted nfw model; weakac — nfw model with weak halo contraction following lAbadi et alj 1 I2OICII) ; 
freeac — nfw model wit h free halo response ranging from adiabatic contraction to expansion; atlas — halo profile adopted by ATLAS3D 
llCappellari et al .l2012al) which is a generalised NFW with fixed scale radius; free — generalised NFW halo with free inner slope (repeated 
from Fig.|4j. Upper: a for bulges (red, solid lines) and disks (blue, long dashed lines); Lower: a for total (bulge+disk) stellar mass. The 
bold histograms show the joint constraints on the mean a. For all models except "ac" the bulges are consistent with a Salpeter IMF. 
For "ac" the bulges are consistent with a Chabrier IMF, but allowing the halo response to vary "freeac" generally finds better fits with 
weaker halo contraction or expansion. 



haloes ( "core" ) , to contracted NFW haloes ( "ac" , "weakac" , 
and "freeac"). The most generic m odel adopted by the AT - 
LAS3D team in their recent study (|Cappellari et al1l2012al ) 
is also included to facilitate comparisons between our works 
("atlas"). 

The contracted NFW halos are meant to represent the 



response of the dark matter halo to gala xy formation. The 
simple model "ac" refers to the standard iBlumenthal et al] 
1 19861 ) formalism, the "weakac" model is the reduced con- 
traction model of lAbadi et al.] (|2010l ). while the "freeac" 
model is a generalisation. In the standard formalism the 
adiabatic invariant is rM{r), where r is galactic radius and 
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M(r) is the spherically enclosed mass within r. Thus 
rf/n = Mi{ri)/Mf{rt), 



(7) 



where ri and rf are the initial and final radii, respectively. 

In order to explore the possi bility of weaker halo 
contraction (e.g., lAbadi et al.l I2OIOI I and even expansion, 
we also consi d er the generalised contraction formula from 
iDutton et all (|2007l ). If the standard contraction ratio 
from Eq. [7] is F = (rf/n), then the modified contrac- 
tion ratio is given by F"^. Standard adiabatic contrac- 
tion llBlumenthal et al.1 [l986l ') corresponds to = 1, the 
iGnedin et al. (12004 ) mode l can be approximated with v 



0.8. the lAbadi et al.1 (|2010l ) model can be approximated with 
ly ~ 0.4, no contraction corresponds to v = 0, while expan- 
sion corresponds to < 0. 

Table [3]lists the dark halo priors for the eight models we 
consider here. We note that the models "core" and "nfw" are 
specific realisations of the model "free" , and the models "ac" 
and "nfw" , are specific realisation of the model "freeac" . 

The posterior probability distribution functions for the 
IMF mismatch parameters of the bulge, disk, and total stel- 
lar mass for all models (our reference "free" model and alter- 
natives) are given in Fig. [7] As expected, the disk masses are 
somewhat sensitive to the choice of dark matter halo. How- 
ever, the bulge masses are remarkably insensitive to the dark 
matter halo model. All of the models, with the exception of 
adiabatically contracted NFW haloes (model "ac"), favour 
Salpeter-type IMFs for the bulges. The "ac" model is the 
only one that favours a Chabrier-type IMFs for the bulges. 
However, the fact that the "freeac" model contains the "ac" 
model allows us to perform a clean model selection proce- 
dure and quantify whether the models with Salpeter IMF 
or the "ac" model provide a better overall description of the 
data. Qualitatively, as can be seen from Fig. |6l the "freeac" 
model prefers values of the contraction index that disfavour 
the standard adiabatic contraction models. Quantitatively 
the comparison between the "freea c" model and the "ac " 
model is given by the evidence ratio (jSivia fc Skillinej[2006l ) . 
The evidence of the "freeac" model is 35 times larger than 
that of the "ac" model", which corresponds to strong evi- 
dence that the fir st model is to be pre ferred according to 
standard criteria (|Sivia fc Skillinell2006l ). Similar consider- 
ations can be made for the IMF mismatch parameter a of 
the disks, although in general it is less well constrained than 
that of the bulge. All models with dark matter prefer disks 
that are Chabrier or lighter, except for the "nfw" and "core" 
models. Those are subsets of the "free" model and are dis- 
favoured by the evidence ratios. 

We thus conclude that it is possible to reconcile 
the data with a Chabrier IMF for bulges if one asserts 
that standard adiabatic contraction is the way dark mat- 
ter halos respond to galaxy formation. However, cosmo- 
logical simulations of galaxy formation often find much 
weaker contracti on than predicted by th e adiabatic contrac- 
tion formalism |johansson et aL I I2OO9I : lAbadi et al.l I2OI0I: 



Tissera et al. 201C ), or even expansion ( Governato et al.l 



2010l : lMaccio et al.H2012D and thus we do not actually expect 
adiabatic contraction to occur in nature. Furthermore, under 
the more general assumption of ignorance about the effects 
of galaxy formation on the underlying dark matter halo, the 
data clearly prefer uncontracted or mildly expanded halos 
with a Salpeter IMF. This is true for a broad range of mod- 
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Figure 8. Effect of external convergence on the posterior distri- 
butions for the IMF mismatch parameter relative to a Chabrier 
IMF, a, for the "free" model. The model "free-|-Kext" (right pan- 
els) has the same dark matter halo as "free" (left panels, repeated 
from Fig. |4]l, except we have included 5% external convergence. 
The effect of external convergence is to lower o by a small amount. 



els, including our "free" models, our "freeac" models, and 
the "atlas" models. 



5.2 Line of sight effects 

Another potential concern is whether excess mass along the 
line of sight could lead us to overestimate the mass associ- 
ated with the main deflector and thus overestimate a. This is 
usually describe d in the lensing literature as extern al conver- 
gence (Koxt, e.g.. lTreu et~al]|2010l : ISuvu et al.ll2010l '). The ex- 
cess/deficit mass along the line of sight acts to first approx- 
imation as a sheet of mass at the redshift of the deflector, 
expressed in units of the critical density Kcrit , and cannot be 
measured with pure len sing arguments due to the well know n 
mass-sheet degeneracy (|Falco. Gorenstein. fc Shapirdll985l ) . 
External convergence affects the lensing observable used in 
our analysis - the amount of mass within the Einstein Ra- 
dius - in a very simple manner. The true mass will just be 
the observed mass obtained by assuming the line of sight has 
the average density of the universe, multiplied by (1-Koxt), 
for small values of Koxt. 

At the relatively low redshift of the SWELLS and 
SLAGS samples, external convergence is very small. Detailed 
analysis of the SLAGS sample shows that the external con- 
vergence is typically a few percent (e.g.. iTreu et al. I l2009l : 
iGuimaraes fc Sodr4l20Tll : ISonnenfeld et al.ll2012i r Thus, we 
expect that external convergence will change our results by 
a few percent at most. For completeness, we repeated all 
our inference assuming Koxt = 0.05. As expected, a shifts 
down by a negligible amount (compare model "free" to 
"free-|-Kext" in Fig. [8]) demonstrating that our inferences are 
robust with respect to line of sight effects. 
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5.3 Cold Gas 

In principle the stellar masses derived from our lensing plus 
dynamics fits are upper limits, because we do not include 
cold gas in our mass models. For the stellar masses of our 
galaxies we expect cold gas fractions of ~ 20 ± 10% (e.g., 
iDutton et allbOllal 'l. roughly equally split between atomic 
and molecular gas. Molecular gas generally traces the stars, 
and is typically le ss than ~ 10% of the s tellar mass, assuming 
a Chabrier IMF (jSaintonge et al.ll201ll ). so the effect on our 
derived bulge stellar masses is expected to be less than 0.05 
dex. The atomic gas typically has a larger scale length than 
the stellar disk, so we expect the atomic gas to subtract mass 
from the stellar disk and dark matter halo. In summary, we 
expect that including observations of cold gas in our mass 
models will not significantly reduce the derived stellar mass 
of the bulges. 



6 SUMMARY 

We have presented mass models of 5 massive gal axies se- 
lected from the SWELLS survey (|Treu et al.ll201ll ') to have 
bulges and disks of comparable stellar mass, as well as star 
forming disks. We combined masses from strong lensing with 
ionised gas kinematics at ~ 1 — 3 disk scale lengths to con- 
strain the parameters of three component mass models con- 
sisting of a bulge, a disk and a generic dark matter halo. 
Our main results can be summarised as follows: 

• The stellar masses of the bulges are well constrained 
by the lensing and kinematic data, independent of the inner 
density slope of the dark matter halo. 

• The bulge masses inferred from the lensing and dynam- 
ical models, Mld , are inconsistent with those obtained from 
the colours, Msps, assuming a Chabrier IMF, but in good 
agreement with those based on a Salpeter IMF. The average 
normalisation of the IMF of the bulges relative to that based 
on a Chabrier IMF, is given by the IMF mismatch parameter 

logio abulge = logio(MLD,bulge/Ms°ps" bulge) = 0.29 ± 0.05 

• The disk meisses inferred from the lensing and dynam- 
ical models are only weakly constrained, due to degenera- 
cies with the dark matter halo, but are consistent with a 
Chabrier-like IMF. The average IMF mismatch parameter 
is found to be logj^Q Odisk = —0.01 ± 0.12. 

• Disks are sub-maximal at 2.2 disk s cale lengths (in 
agree ment with the Disk Mass project, iBershadv et al.l 
l201lf ). However, baryons dominate the potential inside 2.2 
disk scale lengths due to the strong bulge components. And 
thus sub-maximal disks do not imply galaxies are dark mat- 
ter dominated inside 2.2 disk scale lengths. 

• The data marginally disfavour an inner slope of the dark 
matter halo 7 > 1 that would be expected for NFW halos 
contracted according to standard adiabatic contraction pre- 
scriptions. Equivalently, the data favour an uncontracted or 
marginally expanded NFW halo. 

Our main new result is that IMF of bulges of spirals 
is "heavier" than Chabrier, and consistent with a Salpeter 
IMF. Since our data do not strongly constrain the IMF 
of the disks, there are two possible implications. Either 
the IMF varies between spiral galaxies (e.g., massive spi- 
ral galaxies have heavier IMFs than the Milky Way), or the 



IMF varies within spiral galaxies (e.g., bulges have heav- 
ier IMFs than disks). Since all previous constraints on the 
masses of galactic disks seem to rule out Salpe t er IMFs (e.g.. 
Bell fc de J on3 '2OOI'; 'Bcrshady ct al.' '201l|: iDutton et ahl 

Li) 



2011al : Martinssoiii ,2011 : Barna bc ct al.. ,201^ ') we favor the 



latter hypothesis. Even though this result might seem sur- 
prising, it is well known that the stellar populations of bulges 
and disks differ in age and chemical composition. Thus, if 
the IMF reflects the physical conditions at the time of for- 
mation of the stellar populations, it is entirely possible that 
it could be different for the bulge and disk. Clearly, even 
though this is possible, the underlying physical reasons are 
at present unclear. Our hope is that this new piece of the 
puzzle will be a valuable clue for deciphering the mystery of 
the IMF and its variations across the universe. 

To conclude we note that our results are consistent with 
previous work based on completely different techniques and 
samples when re framed in terms of a global IMF (as sum - 
marized by, e.g., iTreu eral]|2010l : rCappellari et al.ll2012bl ). 
Furthermore our hypothesis of a different IMF for bulge and 
disk is consistent with the upper limits on total mass within 
the Einstein radius for the entire SWELLS and SLACS 
sample s, presented in a companion paper (|Brewer et al.l 
l2012bl . in preparation). Whereas SWELLS-III showed that 
the SWELLS and SLACS data are consistent with an IMF 
that changes as a fun ction of galaxy stellar mass of stel- 
lar velocity dispersion, iBrewer et al.1 (|2012bl ) shows that a 
scenario where the IMF is Salpeter-like in the bulge and 
Chabrier-like in the disk is perfectly consistent with the lens- 
ing constraints. 
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